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We consider a four-terminal setup of a two-dimensional topological insulator (quantum spin Hall 
insulator) with local tunneling between the upper and lower edges. The edge modes are modeled as 
helical Luttinger liquids and the electron-electron interactions are taken into account exactly. Using 
perturbation theory in the tunneling, we derive the cumulant generating function for the inter-edge 
current. We show that different possible transport channels give rise to different signatures in the 
current noise and current cross-correlations, which could be exploited in experiments to elucidate 
the interplay between electron-electron interactions and the helical nature of the edge states. 
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Topological insulators are characterized by having a 
nonzero band gap in the bulk while gapless edge modes 
exist on the surfaces pQ [5]. The edge modes are pro- 
tected by time-reversal symmetry and are robust against 
perturbations such as impurity scattering or electron- 
electron interactions. Two-dimensional topological insu- 
lators have been realized using HgTe quantum wells [3] 
and have been shown to display the behavior predicted 
for quantum spin Hall insulators [1]: a strong spin-orbit 
coupling locks the spatial motion of the electron to its 
spin and thus entails the formation of helical edge modes, 
where electrons with opposite spin orientations on any 
given edge propagate into opposite directions. 

In the presence of electron-electron interactions, the 
edge modes of two-dimensional topological insulators can 
be modeled as helical Luttinger hquids (LLs) [SHZ]- This 
extension of the conventional LL theory provides a de- 
scription of the low-energy degrees of freedom and makes 
an exact treatment of the interaction possible. In partic- 
ular, it allows to analyze the nonequilibrium transport 
properties of helical edge modes [5] . Some of these trans- 
port properties have already been investigated in exper- 
iments [iiii]. 

To date, most attention has been devoted to the in- 
vestigation of the average current through edge modes 
as a function of an applied bias voltage [SJIS]. However, 
it has been known for a long time that a measurement 
of the current noise yields additional information, e.g., 
about the nature of the charge carriers, which is not 
accessible from current measurements [TU]. In this pa- 
per, we shall investigate the transport properties of edge 
states of two-dimensional topological insulators using full 
counting statistics (FCS), which yields insights into the 
average current, the current noise, and arbitrary higher- 
order correlation functions of the current. 

We consider the two-dimensional topological insula- 
tor depicted in Fig. [T] The system is insulating in the 
bulk, but the upper and lower edges of the sample sup- 
port counter-propagating gapless helical modes. The 
upper edge contains a right-moving spin-up mode (de- 
noted by the fermionic operator ■0^^) and a left-moving 
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FIG. 1. (Color online) The setup under investigation con- 
sists of a two-dimensional topological insulator with two he- 
lical edge modes. The four corners are connected to electron 
reservoirs held at bias voltages ^jla,b,c,d- Using top gates, the 
upper and lower channels can be coupled locally, leading to 
inter-edge tunneling at a; = 0. 



spin-down mode while the lower edge hosts the 

modes ij^Ri and ?/'Lt- The two modes on each edge are 
coupled by electron-electron interactions, but in a time- 
reversal invariant system, these interactions cannot lead 
to backscattering. The helical LL theory is based on the 
bosonization of the fermionic operators (for a = R,L and 

i^U^) = ^e-^^°'(-), (1) 



where L^qo- denote Klein factors which ensure the correct 
anticommutation relations between the fermionic fields 
[H] and a is a short-distance cutoff. The helical LL the- 
ory can be formulated in different bosonic bases: the heli- 
cal edge basis [13] is the canonical choice and respects the 
separation between the upper and lower edges by using 
the fields (j>i oc (pn^ + ip^^ and (/)2 oc ipj^^ + ipR^. How- 
ever, once tunneling between the two edges is considered, 
the spin-charge basis [13] becomes more convenient. The 
latter mixes states on the two edges and describes the 
system in terms of spin and charge degrees of freedom by 
introducing the spin modes (0s, Og) and charge modes 
{<j)c^ Oc)- For a = R,L = and a =f ,1= -|-, — these 
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are defined by [Tl] 
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(q0c - 6'c + acr(/)s - aOs) ■ 



(2) 



The fields (/>i/ and 9jj [v — c, s) are canonically conju- 
gate, [4>^(x),dyO^{y)] — iTrSfj,i,S{x — y). In terms of these 
fields, the edge state Hamiltonian becomes a sum of in- 
dependent spin and charge terms. The effect of electron- 
electron interactions manifests itself in the Luttinger pa- 
rameters Kc and Kg for the charge and spin sectors. It 
was found that in a helical LL, these parameters are re- 
lated by [U 
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where g = I corresponds to noninteracting fermions and 
< g < 1 for repulsive interactions. 

We assume that the four corners of the sample are con- 
nected to electron reservoirs held at chemical potentials 
IJ'A,B,c,D as shown in Fig. [l] In order to incorporate 
these voltage sources into our theory, we use the ^(a::)- 
model |15H17j . In this approach, the interacting region 
of the ID wires extends from x — —L/2 to x — L/2 where 
L is the sample length. The electron reservoirs, on the 
other hand, are modeled as noninteracting wires located 
at I a; I > L/2. This approach has been used successfully 
to calculate current and noise properties of LLs [18] . In 
the spin-charge basis, the Hamiltonian describing the in- 
teracting edge states thus becomes Hq — Hc + Hs, where 
(using e — h — 1 throughout this paper) 
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where g{x) = g for \x\ < L/2 and g{x) — 1 for \x\ > L/2. 
The voltage sources shift the chemical potentials in the 
noninteracting parts of the wire and affect both spin ori- 
entations in the same way. For instance, the voltage fiA 
leads to a term fj,A dx[pii]-{x) + pL^{x)]. Analo- 

gous expressions proportional to pb.c,d appear at the 
three remaining corners. Using bosonization, the voltage 
sources can be described by the Hamiltonian 
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dx[E,{x)^,{x)+Es{x)e,{x)]. (6) 



Note that the applied voltages couple to the charge field 
<f>c but to the conjugate spin field Og- This reflects the 
nonequilibrium spin-charge duality found in Ref . [5j . The 
voltage terms in the charge and spin sector are given by 

E,^,{x) = (mc ± iiB)5{x - L/2) - [pD ± Pa)5{x + L/2.). 

(7) 



Finally, we assume that near the point a: = 0, the setup 
is constricted in such a way that the helical modes on 
both edges spatially approach each other. In this case, 
local tunneling becomes possible. The possible tunnel- 
ing terms in the presence of interactions in the wire 
were identified in Ref. [7]. First, spin-conserving single- 
particle tunneling leads to 



(8) 



where V'qo- = V'ao- (x = 0) . In addition, there may be tun- 
neling of either charged or spinful particle pairs, which 
gives rise to the following terms 



(9) 
(10) 



An analysis of the scaling dimensions of these three terms 
reveals that is the dominant contribution for weak in- 
teractions (5 ~ 1). Towards stronger interactions, how- 
ever, their magnitudes may becomes comparable and it 
is generally no longer justified to neglect pair tunneling 
[7]. For 1/2 < 5 < 2, all terms remain irrelevant in 
the renormalization group (RG) sense such that pertur- 
bation theory in T^^c.s is applicable. Note that using a 
mapping between the limiting cases of weak and strong 
inter-edge tunneling, our results can also be applied to 
the case where the sample is almost pinched off and only 
residual tunneling connects the left and right sides [7|. 

Full counting statistics (FCS) is a convenient tool to 
extract information about the nonequilibrium transport 
properties of this system. The objective of FCS is the 
calculation of the cumulant generating function (CGF) 
Inx(A), which in its simplest form yields the average 
current, the zero- frequency noise and other current cor- 
relation functions. We are interested in calculating the 
statistics of the tunnel current, so the CGF is defined as 



Inx(A) = In ( exp 



(11) 



where 5Naa denote the number of particles which have 
been scattered out of the edge mode V'ao- during the mea- 
surement time T. For large T, the CGF becomes pro- 



portional to T. Derivatives of Eq. (Ill with respect to 
the counting fields \aa then allow a calculation of the 
average tunnel currents (SLaa) = (SNaa-) /T and of arbi- 
trary correlation functions, e.g., {{SNarrSNa'a')) which is 
proportional to the zero-frequency current noise. 

A direct calculation of the CGF using the definition 
(11) is difficult. It helps to take advantage of its rep- 



resentation as a time-ordered expectation value on the 
Keldysh contour [T5], 



Inx(A) = In ( Tc cxp 



i J dsT^{s) 



(12) 
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Here, Tq denotes the time-ordering operator on the 
Keldysh contour C, which consists of two branches: C_ 
reaches from — oo to +00 whereas C+ leads back to —00 
PO] , The subscript V in the expectation value indicates 
that it is taken with respect to the ground state of the 
Hamiltonian Ho + Hy- In order to construct the tun- 
neling operator T^, all fermionic operators have to be 
furnished with counting fields, 



(13) 



The counting fields have a nontrivial time-dependence: 
Ka{t) = ±A„a for t e C± [H]. Then, can be con- 



structed by using the substitution ( 13 ) in the definition 
of the tunneling Hamiltonian T^+Tc+Ts, and bosonizing 
it using Eq. ([T]) . 

Unfortunately, an exact calculation of the expectation 
value in Eq. \12\ in nonequilibrium and in the presence of 
electron-electron interactions is not possible. Therefore, 
we consider the tunneling as a small perturbation. Then, 



we can use a linked-cluster expansion of Eq. ( 12 ) and 



truncate the series after the second (leading) order in 
7e,c,s- Up to this order, the total CGF remains a sum of 
three contributions due to the individual tunneling terms. 



In X = In Xe + In Xc + In Xs 



(14) 



The single-particle term as well as the terms due to tun- 
neling of charged pairs and spinful pairs are given by 
(using V = c,s) 



InXe(A) 



llT 
(27ra)2 



' m=± a=^,i=± 



X Cc,^{ms)Cs,ma{ms), (15) 

2q- poo 

Inx.(A) = / ds ^ (e™^" - l)C,.M^s), 



and depend on the following linear combinations of the 
counting fields, (ct =t,|= +, — ) 



(16) 



The CGF is determined by the bosonic correlation func- 
tions Ciy^(s). For v — c,s and C G these are defined 
by 

a,c(i) = ^e^^-^'^We-^^^-^"))^ . (17) 

The calculation of these expectation values is a nontriv- 
ial task due to the position-dependence of the Luttinger 
parameter g{x) and due to the presence of the voltage 
sources. However, the necessary techniques for ordinary 
and helical LLs have been developed and are explained 
in detail in Refs. [51 [TH]. The correlation functions are in 
principle known for arbitrary system length L and even 



at nonzero temperatures. Finite-length effects in the tun- 
nel current have been investigated in Ref . [8j . They have 
been found to lead to an oscillatory current-voltage char- 
acteristic which is observable only at very low bias volt- 
ages. Therefore, in order to keep the results compact, we 
shall focus on the case of zero temperature and infinite 
length. In this case, the correlation functions become (for 
v = c.s) 



Cvxit) = exp 



it) 



(18) 



where Uc ~ vp/a is a high-energy cutoff and = 
l/Kg — g. The effective voltages fic^s affecting the charge 
and spin sectors are given by 

fJ-cs = l:{fJ.c T fJ-D ± fJ.B ~ fJ'A)- (19) 



The time-integration in the CGF ( 15 ) can be performed 



exactly and one finds the main result of this article: up 
to the second order in tunnel amplitudes, the CGF of the 
tunnel current consists of the terms. 



i^xe(A) = _„,,,^ 7:;^, ^ 



7ra2r(.g + 1/g) 



{flc + (^fJ-s) 
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(20) 



Hence, the CGF is a sum of three statistically indepen- 
dent Poissonian terms. The fact that all contributions are 
Poissonian is reasonable: by using perturbation theory in 
the tunneling, we assumed that inter-edge tunneling is a 
rare event. Next, we will discuss physical predictions for 
the current and current correlations which can be ob- 
tained from Eqs. (Il4| and (pO). 



The spin-resolved tunnel current can be calculated 
using {5h) = {dN,)/T = (l/r)51nx(A)/9(zA,)|A=o. 
Most notably, the three transport processes lead to differ- 
ent power-law current-voltage characteristics whose ex- 
ponents are different combinations of the Luttinger pa- 
rameter g. Power laws are a common feature of trans- 
port in interacting systems [5T]. The single-particle cur- 
rent has already been found in Ref. [8^. For small in- 
teractions g ~ I, single-particle tunneling becomes lin- 
ear in the applied voltage, and is thus dominant because 
pair-tunneling is cubic in /j,c,s- Charged-pair tunneling, 
with an exponent Ag—l, becomes more relevant towards 
stronger repulsive interactions (1/2 < g < 1). The ex- 
ponent for spinful pair tunneling, on the other hand, is 
4/(7 — 1, so its contribution is smaller than the other pro- 
cesses for repulsive interactions. All power laws are con- 
sistent with an RG analysis of the tunneling terms [7]. 
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The fact that ah three contributions to the average cur- 
rent depend on different combinations of the apphed voh- 
ages iJLc,s could aUow an experimental distinction of the 
different transport processes. 

The noise of the tunnel current is proportional to 
the second moment of the CGF, = {{SN^SN^)) = 
{l/T)dHnx{X)/d{iX^)^\x^o- It is given by 



na^Tig + 1/g) 



(27r)3a4r(4g) 
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(27r)3a4r(4/.g) 



(21) 

4/9-1 
4/ff 



Similar to the average current, the contributions due 
to single-particle tunneling, charged pair tunneling and 
spinful pair tunneling couple differently to the applied 
voltages fj,c and fis- In the absence of scattering and at 
low temperatures, the transport in clean ID LLs does 
not produce zero-frequency shot noise 18]. Therefore, 
the tunnel current noise can directly be measured as the 
noise in either of the four contacts. 

An even clearer distinction between the transport pro- 
cesses can be obtained by considering cross-correlations 
of the currents along different edges. These can be re- 
garded as fermionic Hanbury Brown and Twiss (HBT) 
correlations and were shown previously to be significantly 
modified by interactions in multi-terminal setups |22j . 
Such correlations have already been measured experi- 
mentally in quantum Hall bars [23] . In the present setup, 
the interplay between interactions and the helical struc- 
ture of the edge modes makes the result particularly in- 
teresting. Let us illustrate this with a simple bias con- 
figuration, — ^2y, whereas fJ,B,c,D = 0, such that 
i^-c — fJ-s = V- We find 



5hbt = {{SN^SN^)) 
1 / 



(22) 



(27r)3a4 \r(4/5) 



U!c 



r(45) 



1 ^ 



The single-particle tunneling does not show up in the 
HBT correlation because it does not couple spin-up and 
spin-down electrons. In contrast, particle-pair tunneling 
does yield HBT correlations. The sign of the correla- 
tions is positive for charged-pair transport and negative 
for spinful-pair transport. The physical reason for this 
behavior can be understood from the different forms of 
the tunneling terms, see Fig. [2j 

Depending on the interaction strength, either charge 
transport or spin transport will be dominant [7|. In 
the noninteracting limit {g — 1), S'hbt is proportional 
to |T^|^(7s ~ 7c) and thus allows a direct comparison 
between the strengths of two pair-tunneling processes. 
For repulsive interactions, 1/2 < g < 1, the different 
power laws associated with both processes will lead to 
a nonmonotonic S'hbt (1^)- These features of the cur- 
rent cross-correlation are therefore direct evidence of pair 
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FIG. 2. (Color online) Graphical depiction of charged-pair 
tunneling (Tc) and spinful-pair tunneling (T3) according to 
Eqs. ([9|-( |10[ |. Due to the helicity, a single Tc event turns, 
e.g., two right-movers into two left- movers and thus changes 
the charge transferred from left to right by two. Ta affects 
the spin sector in an analogous way. Both processes give rise 
to cross-correlations between spin-up and spin-down tunnel 
currents {{SNi^SNi)) . As Tc transports two particles in op- 
posite directions, whereas Ts transports two particles in the 
same direction, the sign of the cross-correlations is different 
for both processes. 



tunneling and can be used to distinguish different trans- 
port regimes experimentally. Due to the helical struc- 
ture of the edge modes, an experimental measurement of 
S'hbt at this bias configuration does not require a spin- 
resolved current measurement. Instead, it can be deter- 
mined by measuring the cross-correlation of the counter- 
propagating spin-up and spin-down currents on the lower 
edge, Shbt oc {{iLtlRi)) ■ 

In realistic systems, spin-orbit coupling may also cause 
spin- flip tunneling processes [2l[8], i.e., tunneling terms 
of the form An analysis of such processes us- 

ing FCS reveals that they do not cause correlations of 
HBT type. The physical reason is that spin-flip tunneling 
represents a single-particle process, whereas according to 



the definition (22 1 correlated tunneling of two particles 
is required to produce HBT correlations. Therefore, the 
results for S'hbt remain valid even in the presence of spin- 
flip tunneling. Moreover, all features should be visible for 
weak interactions, where the amplitudes for charged-pair 
and spinful-pair tunneling are comparable. 

In conclusion, we have investigated transport in heli- 
cal edge states of a two-dimensional topological insulator. 
For weak local tunneling between the upper and lower 
edges, we used perturbation theory to calculate the cu- 
mulant generating function for the tunnel current. This 
allowed us to analyze the current noise and current cross- 
correlations which emerge in this setup. Different trans- 
port processes produce strikingly different terms in the 
cross-correlations. This could be exploited to experimen- 
tally distinguish different transport regimes and investi- 
gate the interplay between electron-electron interactions 
and the helicity of the edge modes. 

The author wishes to thank A. Komnik for valuable 
discussions and acknowledges support by the Swiss NSF. 
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